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Abstract Based on Newton's laws reformulated in the Hamiltonian dy- 
namics combined with statistical mechanics, we formulate a statistical me- 
chanical theory supporting the hypothesis of a closed oscillating universe. 
We find that the behaviour of the universe as a whole can be represented by 
a free entropic oscillator whose lifespan is nonhomogeneous, thus implying 
that time is shorter or longer according to the state of the universe itself 
given through its entropy. We conclude that time reduces to the entropy 
production of the universe and that a nonzero entropy production means 
that local fluctuations could exist giving rise to the appearance of masses 
and to the curvature of the space. 
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1 Introduction 

Cosmology, as an important discipline of scientific knowledge, tries to ad- 
dress fundamental questions beyond the realm of sciences and influencing 
our philosophical and even religious comprehension of the world. Perhaps 
the major question posed by cosmology is the origin of our universe. Nowa- 
days, the most important theory providing an answer to this question is 
the Big-Bang theory [Tj. Aside from others details, we are mainly interested 
here in one of the implications of the Big-Bang theory, i. e., "that the uni- 
verse will one day end" , [T] . According to a current speculation, in a closed 
universe, after the expansion eventually stops, a contraction will follow lead- 
ing to an implosion into a singularity a process known as the Big Crunch. 
This oscillating behavior of the universe has important implications over 
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the concept of time, first because it is not clear that the initial and final 
states of the universe can be correctly conceptualized through the notion 
of time, understood as a fundamental quantity, and secondly because this 
behavior suggests time is created and destroyed at the beginning and the 
end of the universe. 

In contrast with these implications, we know that in all fundamental 
physical theories time plays a central role. Mechanics, quantum mechanics 
and relativistic theories require time as a quantity through which the evo- 
lution of the system is conceptualized. However, thermodynamics does not 
require the concept of time in order to establish the "direction" in which 
spontaneous natural processes occur, that is, towards the state of thermo- 
dynamic equilibrium. These considerations lead to the following question: 
Is the time a fundamental quantity to describe the evolution of a system 
during a spontaneous natural process? 

Beside the Big Bang theory other alternatives with respect to the origin 
and nature of the universe have been posed: is it open or closed...? it is 
cyclic?, or is there an arrow of time related with entropy and irreversibility?. 
Therefore these questions are still a recurrent topic of debate [2]-|10j. 

Here, we present a cosmological statistical mechanics theory of a closed 
universe not in equilibrium, oscillating between two thermodynamic equi- 
librium states which runs parallel to other cosmological models proposing 
an oscillating or cyclic universe [3] , [1] , [TO] ■ Our statistical mechanical the- 
ory might complement these models. We assume that physical interactions 
give rise to time and space, that is, we consider that if physical interactions 
did not exist, time and space would not exist either. Our theory provides 
us with an expression of the entropy production of the universe based on 
first principles from wich one infers that time is inhomogeneous and re- 
duces to the production of entropy. Moreover, our theory implies that there 
is no global time arrow although it explains the existence of irreversible lo- 
cal phenomena [llj as the corresponding ones in a nonequilibrium fluid. A 
more general statistical theory can be elaborated by incorporating quantum 
and relativistic effects through a description in terms of the set of reduced 
density operators instead of the distribution vector in addition to the gen- 
eralized von Neumann equation and the generalized von Neumann entropy. 
Density operators will account for quantum and also relativistic effects since 
the metrics of the space-time will also be contained in this description [12]. 
However, since the stochastic dynamic is determined by the eigenvalues of 
the Hamiltonian, this more mathematically elaborate theory will essentially 
lead to the same conclusions. 

To mathematically formalize the previous statements, we will consider 
an isolated system composed by N particles or bodies (the closed universe) 
whose large-scale structure is governed by long-range interactions falling 
into the Hamiltonian of the universe itself. 

In classical mechanics the state of an N-body system at any time is 

given by a set of 3N generalized coordinates , q%N and 3N conjugated 

generalized momenta pi, ,P3N- The value of these 6N variables defines 
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a point in the 6N-dimensional phase space T corresponding to the system. 
This representative point of the system moves in the phase space along a 
trajectory determined by Hamilton's equations 

dH dH 

dpi dqi 

where H is the Hamiltonian of the universe. Therefore, Hamilton's equations 
describe a flow in T space whose density flow F(T, t) varies according to 
the Liouville equation 

d 



i-F=LF (2) 



with 



\ dqi dpi dpi dqi J 



(3) 

• \ 1 mi 1 n)i t tin 1 it\i 1 

being the Liouville operator. The description in terms of the Liouville equa- 
tion is completely equivalent to that in terms of Hamilton's equations |13j . 
If the system occupies a volume V in the configurational space the density 
will vanish out of this volume, a fact which can be used to show that the 
Liouville operator is Hermitian [14] . 

Given an initial density F(T,0), the formal solution of the Liouville 
equation ([2]) can be written 

F(T,t) = e~ ltL F{T,0). (4) 

Since L is Hermitian, all its eigenvalues are real [15], which according to Eq. 
(j4|) implies that F(T, t) will have an oscillatory behavior. This is an impor- 
tant general result based on the fully microscopic Hamilton dynamics of the 
system which allows a rigorous formulation a closed universe in terms of 
the Liouville equation A consequence of the Hermitian character of the 
Liouville operator is that one may formulate a statistical mechanical theory 
of a closed oscillating universe in which the evolution is parametrized by 
means of the entropy of the universe itself, and in turn time is parametrized 
through the rate of variation of the entropy, i. e., the entropy production. 

The paper is organized as follows. In section 2, we formulate the Gen- 
eralized Liouville equation and give its formal solution. In section 3, we 
postulate the nonequilibrium entropy and obtain the entropy production of 
the universe itself. Section 4 is devoted to the derivation of the equations of 
the entropic oscillator corresponding to the universe and finally in section 
5 we discuss our main conclusions. 



2 Generalized Liouville equation 

Let us consider that the Hamiltonian of the system is given by 
N p 2 1 N 
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with rrij being the mass of a particle, q_, the position vector of the j — th 
particle and Pj its conjugated momentum. Moreover <p (|cjj — qfc|) = 4>j u l& 
the interaction potential. The complete statistical description of the system 
can be given in terms of the distribution vector f [T3] , [TJJ 

f(t) = {f J 1 (x 1 ,t)J 2 (x 2 ,t), ,f N (x N ,t)} (6) 

which is the set of all the n-particle reduced distribution functions, 

fn = J F(x N ,t) dx n+1 ...dx N , (7) 

where Xj — (q^, Pj) and x n = {x\, x n }, n = 0, ,N. Here, f Q — 1 and 

Jn(x n ) = F(x N ,t). The dynamics of the distribution vector f is obtained 
by integrating the Liouville equation |(5J according to the definition of the 
rt-particle reduced distribution functions 0, which gives 



d id 

— f n = {H n , /„} + (N ~ n) / Fj, n+1 — f n+1 dx n+1 . (8) 
Here, F J>+1 = -Vj0 J>+1 and 



n 2 i n 

is the n-particle Hamiltonian. The first term on the right hand side of Eq. 
([8]) represents a Hamiltonian flow, while the second, a non-Hamiltonian con- 
tribution due to the coarse graining of the description. The set of equations 
represented by Eq. ([8|) constitutes a hierarchy of coupled equations, the Bo- 
goliubov, Born, Green, Kirkwood, Yvon (BBGKY) hierarchy of equations 
[T5] . [TT] representing the full microscopic description of the system. The 
BBGKY hierarchy, which can be represented in a compact way 

i^f(t) =£f(t), (10) 

constituting the Generalized Liouville equation. By comparison of Eqs. (|8| 
and (fTTJ|) it can be seen that the Generalized Liouville operator C splits into 
a diagonal Hermitian part VC defined through [T3], [TT] 

(n |P£| n') = i [H n ,f n ] P 8 n ,, n , n > , (11) 

where |n) represents the n-particle state, and a nondiagonal non-Hermitian 
part QC 



(n\ QC \n') = i | (N - n) ^ F j>+ i-^/„ + i(It„ + i | 



+i , n > 1. 

(12) 
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Therefore, Eq. (fTOjl can be rewritten 

i^f{t)-VCf{t)=QLf(t) (13) 

which explicitly manifests the Hermitian and non-Hcrmitian contributions 
to the dynamic of f(i). In Eq. (fT3|) , the term VCf(t) plays the same role 
as the Liouville operator in Eq. ([2|) and since they are both Hermitian, 
they induce comparably incompressible flows. On the other hand, the term 
Q£f(t) introduces long range correlations in the dynamics, which are re- 
lated to the effects of coarse-grained interactions, and thus to dissipative 
effects in a reduced description. This implies that Eq. (fl3|) provides us with 
more information than Eq. l[2]). In fact, as we will see, the non-Hermitian 
term is responsible for the approach to equilibrium, constituting a proof 
of the adequate description in terms of f (t) instead of the full phase-space 
distribution function F. 

The formal solution of Eq. (p~3j) is 



f (t) = cxp (-iVCt) f (0) + exp {iVCt) [ dr exp {-iVCr) (-iQ£)f (r) , 

Jo 

(14) 

which, proceeding iteratively, reduces to 

f(t)=U(t,0)f(0) , (15) 

where the evolution operator U (t, 0) is given by a perturbative development 
as 

U(t,0) = y2 / dti / dt 2 / dt 3 / dtjX 

j=o Jo Jo Jo Jo 

V(i,ii) V(t j - 1 ,t j )exp(-i'P£tj) ■ (16) 

Here, V tj) = cxp [iVC (tj-i — tj)] (—iQC) are nonHermitian propa- 
gators, tj < tj-i < < <i < to = t, and the integration proceeds from 

right to left. Since VC is Hermitian, all its eigenvalues are real [TS], which 
means that f (t) as given through Eqs. (| 1 5|) and (|16p will have an oscillatory 
behavior. This oscillatory behavior has strong consequences in the evolution 
of the closed system, as we will show in the next sections. 



3 Nonequilibrium entropy 

In order to establish a connection between the N-particle microscopic de- 
scription given by the generalized Liouville equation (I13p and a macroscopic 
coarse-grained description, it is necessary to consider the statistical defini- 
tion of the entropy. To achieve this objective, we start by noticing that the 
Gibbs entropy 



S N = -fcsTr(FlnF) = -k B / FlnFdx 



(17) 
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where ks is the Boltzmann constant, is a constant of motion under the Li- 
ouville dynamics given through Eq. |2]) , hence alternatively, as the nonequi- 
librium entropy for the N-body system, we propose 

S = -k B Tr{f]nf- q 1 f}+S eq 

= / fn^-r^— dx i dx n + S eq , (18) 

n=1 J Jeq,n 

a convex functional of the distribution vector which generalizes the Gibbs 
formula. In Eq. (|18p. S eq is the thermodynamic entropy (i.e. the equilibrium 
entropy) and f eq is the equilibrium distribution vector satisfying £f eq = 0, 
the Yvon-Born-Green (YBG) equilibrium hierarchy [T5], Therefore, f eq is an 
eigenfunction of C with eigenvalue 0. 

From the convexity of the logarithmic function, it has been proven [16] . 
[TT] that the entropy S varies between the bounds determined through 

> S eq > S > S N + S eq . (19) 

Making use of Eq. (|10[) we can obtain the rate of change of S, which is the 
entropy production 

dS 



. u - ? feTr{£fln(f- 1 f)} 



yEE / to - k B T gz- + Yl F ^ + ( N - «) ^ dxn 

n=i i=i J \ i-t-i— 1 / 



n=l j=l \ J 

(20) 

Here, Tj(x n ) is the force on the j-th particle from the N — n particles not 
contained in the cluster of size n, and is defined through the relation: /„ (x n ) 
Tj(x n ) = J Fj .„ + i/„ + ida; n+ i, and T is the kinetic temperature obtained 
by taking into account that the dependence of f e q,n on the velocities is 
given through a local Maxwellian. The entropy production given in Eq. (|20p 
vanishes at equilibrium and in any other case it should not be necessarily 
zero. In addition, because Pj is arbitrary, 

53 F hl + {N -n)^ = k B T— ln/ e ,,„ (21) 

is sufficient to satisfy the extremum condition SS/Sf n \ eq = 0, with S = 
dS/dt. Precisely, Eq. (|2"Tj) gives rise to the YBG hierarchy previously men- 
tioned [TT],[T7], and expresses a balance of forces: in the right-hand side of 
this equation there appear the mean force [17] and in the left-hand side the 
sum of two terms; the first is the force due to the van der Waals interac- 
tions with the n — 1 fixed particles different from the j-th particle in the 
n-th cluster (responsible for the compressions and dilatations of the n-th 
cluster), while the second which introduces long range correlations is the 
average force on the j-th particle from the remaining N — n particles of the 
system. 
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4 Entropic oscillator 

In Section 2, we have demonstrated through Eq. (fTS"|) that the distribution 
vector f(t) corresponding to the closed universe possesses an oscillatory 
behavior and in Section 3, we have established the relation between f(t) 
and the entropy S through Eq. (fT8|) . Therefore, as a consequence of these 
S(t) will also have an oscillatory behavior which can be described by 

if one assumes that S(t) behaves as a free oscillator with an "elastic" time 
dependent coefficient K e f f (t) defined through the relation 

Keff{t) = \w S - (23) 

A characteristic value of K e f f (t) can be estimated from the consideration 
that the period of the entropy oscillations of the closed universe should be 
coherent with the Poincare cycles, t p , which leads to 



2tt 



" yjK eff (t)dt. (24) 



From these considerations it follows that the entropy of the closed universe 
will oscillate between two "equilibrium" states: > S eq > S > Sn + S eq . 
This is one of the main results of this paper, since it strongly suggests that 
our statistical mechanical theory of a closed universe might complement 
the standard cosmological model of an oscillating universe based on purely 
general relativity grounds [1] . 

The implications of this theory on the concept of time can be extracted 
by first noticing that the dynamics (f2"2"|) for the entropy of the universe 
corresponds the Hamiltonian [18] 

n(s,S,t)=U ) +#(5), (25) 



2 VV^TTW, 

where $(S) = (l/2)(5-5*) 2 ,with S* = S eq + S N /2. This indicates that the 
central quantities in the statistical description are related to the entropy S 
and its rate of change S. Moreover, since the solution of Eq. ([22| is given 

by 

S(t) - S* = A s cos [z{t) + a] , (26) 

where As = Sjv/2, and a is a parameter containing the initial conditions, 
and where the rescaled time z(t) is defined through the relation 

*(*) = / yjK eff (t>)dt', (27) 
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the solution given through Eq. (|26|) of the Eq. (f22|) enables us to parametrize 
time in terms of the entropy through the relation 

z(S) = ±^= { - dt ' (28) 

For the potential <P(S) already introduced we explicitly obtain z(S) = 
arcsm 



(S — S*)/\f27i . From this relation we may conclude that time 
depends upon the amount of entropy generated by the universe during its 
evolution towards the final equilibrium state. More precisely, in view of Eq. 



([25]) . the quantity y/H — <P(S) constitutes the scaled entropy production 
and therefore time in the universe is self-generated by interactions through 
the entropy production. From this analysis it follows that during the evolu- 
tion of the universe, time is non-homogeneous due to the non-homogeneous 
character of the entropy production of the universe itself. 



5 Conclusions 

In this article, we have proposed a statistical mechanic theory of a closed 
universe configured by physical interactions. These interactions define the 
lifespan of the universe itself and its volume in the configurational space. 
The whole universe can be represented as a point in a two-dimensional 
phase space determined by the entropy S and the entropy production S. 
The universe evolves in its phase space along a closed trajectory which 
constitutes the phase portrait of a free entropic oscillator. Since the entropy 
production depends on the state of the universe through its instantaneous 
entropy, we conclude that time has a nonhomogeneous character. 

The key point in our theory lies in the definition of the state of the 
universe by means of the distribution vector and the postulate of the Gen- 
eralized Gibbs entropy as a functional of the distribution vector. Unlike the 
Gibbs entropy which is a constant of motion under the Liouville dynamics, 
the Generalized Gibbs entropy is a varying quantity whose variation is given 
by the Generalized Liouville dynamics. 

Two main results follow from our analysis: The first one is that our 
theory predicts a closed oscillating universe, in accordance with the full 
microscopic description given by the Generalized Liouville equation, and 
that time is a consequence of dissipative interactions that may be quantified 
by the entropy production of the universe itslef. 

Quantum and relativistic effects can be taken into account by formu- 
lating our theory in terms of the set of reduced density operators and the 
corresponding generalization of the von Neumann equation with the suit- 
able metric, which will be addressed in further works. As it was shown, 
since the main results of the present work are a consequence of the nature 
of the Liouville eigenvalues, despite of the details of these more convoluted 
mathematical formulations, the main conclusion of our model will still be 
valid. 
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